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Changes in the propagation of a VLF signal clue to the changes in phase velocity caused 
by a small localized depression in the height of the ionosphere along or near the line joining 
the transmitter and receiver are investigated using elementary diffraction theory. It is 
shown that both the phase and amplitude of the received signal may be changed if the hori- 
zontal extent of the depression is great enough and if it is sufficiently close to the path. 
If the depression is at greater distances from the path, the amplitude and phase of the 
received signal merely oscillates about the undisturbed values. An account is given of 
the use of Cornu's spiral in estimating the effects of a particular disturbance. 



1. Introduction 

It is sometimes of interest to be able to estimate 
how the phase and amplitude of a VLF radio signal 
will be altered if the height of a small area of the 
ionosphere, along or near the path of the radio wave, 
is altered. Such ionospheric height changes, in the 
form of depressions, might be caused, for example, 
by the localized "dumping" of electrons from the 
radiation belts or by a small nuclear explosion at 
heights such that the resulting prompt gamma and 
x radiation produce ionization at limited ranges. 

Now, in practice, the signal from a distant trans- 
mitter is comprised of energy received over a finite 
range of azimuth angles and not merely that received 
along the line between receiver and transmitter. 
Thus, a disturbance in the ionosphere must subtend 
an angle at the receiver which is comparable with 
the azimuthal sector in which the signal is received 
in order to produce an appreciable effect on the 
received signal. Otherwise energy will be diffracted 
around the disturbance, and it will cause little effect. 

Again, the width of the wavefront which makes 
the significant contribution to the received signal 
depends on its distance from the receiver and 
transmitter and is a maximum when it is equidistant 
from them. Thus, if the ionospheric disturbance is 
sufficiently small and of fixed dimensions, the effects 
it produces will depend on its location along the 
path. It will be seen, therefore, that it is not 
permissible, in general, to apply an earlier result 
obtained by Wait [1961] for an infinitely wide 
disturbance to a disturbance which is of limited 
width (perpendicular to the path) . 

The situation is further complicated, since a dis- 
turbance of the type outlined above, occurring at 
random, is more likely to occur on one side or the 
other of the path rather than on it. Nevertheless, 
in such a case, the propagation of the signal is still 
likely to be affected, although only slightly if the 
disturbance is sufficiently remote from the path. 



In this paper, these effects will be investigated by 
using an approximate method based on Fresnel 
diffraction theory. In order to clarify the work, 
several simplifying assumptions will be made. First, 
it will be supposed that propagation is two dimen- 
sional (in the horizontal plane), and that the earth is 
Hat. It will also be assumed that all distances are 
sufficiently large that only the first order waveguide 
mode need be considered. The effect of a local 
change in height of the ionosphere may then be 
associated with a local change in phase velocity of 
the single propagating mode |Wait, L961]. This 
assumption is reasonable if the slope in the ionosphere 
is small. In any case, it will be assumed that the 
depression does not generate significant higher order 
modes and that refraction can he ignored. Finally, 
it will be assumed that even though the ionospheric 
disturbance is of finite size, its effects can be repre- 
sented by a "lumped" phase shift on that portion 
of the wavefront which passes through the center 
of the disturbance 4 and which is bounded in the 
transverse direction by the edges of the disturbance. 

The value of the phase shift introduced by a par- 
ticular change in the ionosphere depends of course on 
the change in phase velocity which is produced. 
This then depends on the size and shape of the iono- 
spheric disturbance. In general there is an infinite 
combination of sizes and shapes of disturbances 
which could produce the same lumped phase shift. 
The special case of a disturbance in which the height 
of reflection of VLF radio waves cosinusoidally with 
distance has been treated by Wait [1961 1. Thus, in 
this paper it will be assumed that the actual phase 
change <p suffered by a ray of infinitesimal width in 
passing under a particular ionospheric depression 
has already been determined, by using methods such 
as that discussed by Wait [1961 |. 'Phis paper will 
be solely concerned with investigating the effects of 
changes in the width and the position of the dis- 
turbance. In particular the effects of a disturbance 
on the phase and amplitude of the received signal 
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will be calculated for two specific cases: (1) when 
the center of a disturbance of varying width is on 
the great circle path, and (2) when a disturbance of 
fixed width moves across the path. It will also be 
shown that the effects of a small circular disturbance 
tend to be proportional to the square of its diameter. 
Thus small disturbances are much less likely to pro- 
duce observable effects than an inappropriate or 
naive application of Wait's [1961] calculations would 
indicate. Finally, a consequence of including the 
effect of the earth's curvature will be discussed. 

A related problem has been discussed by Wait 
[1961]. He was primarily concerned with the effects 
of a disturbance which, while of finite dimensions 
along the direction of propagation, was infinitely 
large in the transverse direction. Wait was pri- 
marily concerned with the change in phase delay 
which would be expected as the longitudinal length 
and depth of the depression were varied. He found, 
inter alia, that the phase change produced by a de- 
pression of this form was independent of its position 
between the transmitter and receiver. Recently, 
Wait [1964 a, b] has extended these results to cover 
the case of a disturbance of finite width as it moves 
across the path, and he has presented calculations for 
a disturbance of both semi-infinite width and finite 
width. In these later papers, he shows that the 
effects on the received signal depend on the distance 
of the disturbance from the receiver and transmitter. 
This same subject is treated in this paper, using a 
somewhat different approach. 

2. Application of Fresnel Diffraction Theory 

In this paper the approach is that used in text- 
books on optics [e.g., Wood, 1934] in discussing dif- 
fraction by a thin wire or other opaque object. 
Here, however, the analogous diffracting object will 
be transparent rather than opaque, and will intro- 
duce a phase delay which differs from that which 
would be obtained in its absence. 

In this application, it will be assumed that each 
elementary portion, ds, of the initially cylindrical 
wavefront emitted by the transmitter reradiates 
uniformly in all directions. Thus, the signal ob- 
served at the receiver is the sum of the signals 
reradiated by each element of the wavefront. The 
effect of the disturbance is to cause a phase change 
in that portion of the wavefront covered by the dis- 
turbance. The first step is to determine the phase 
of the signal reradiated by an element of the wave- 
front, ds, situated at a distance, s, from the line 
joining the transmitter to the receiver. The geom- 
etry of the situation is shown in figure 1, in which 
T and R are the transmitter and receiver, respec- 
tively. The phase disturbance is centered at P, 
and PQ is a portion of the cylindrical wavefront from 
the transmitter. Thus, TP and TQ are equal to a. 
The distance QR is b, while PR is b + 8. PU is a line 
perpendicular to TR drawn through P. 

The geometry of the figure shows that 

PR 2 =UR 2 +PU 2 




Figure 1. The geometry of the propagation* of a wavefront VW 
from the transmitter T to receiver R. 



or 



Since 



(b+s) 2 =(a J rb— a cosa) 2 +(<z sin a) 2 



a^s/a 



8=s 2 (a+b)/2ab. 
Then the field at R, due to ds. is 



(1) 



E=E sin r^-y 0+6 + 5)1 ds, 

where \E \ is a constant relating the observed field to 
the transmitted power, etc. 

Thus, the total field due to the portion of the wave- 
front between V and W, for which the distances, s, 
from Q are s 2 and s u is given by 



E t =E ( * 2 sin \<»t~ (a+b+8)l ds 



(2) 



which on expanding the sine term can be written in 
the form 



where 



I! 



and 



E t =RE sm[ wt-~ (a+b)-6~\> (3) 

HX cos ^]HX sin ^T (4) 
P 

_i *^ 5 i 



0=tan 



2tt<5 , 
sm — as 

A 



J" 



2tt<5 , 
cos — ds 

A 



The integrals are Fresnel's integrals and can be 
evaluated from tables, as discussed later. 

It should be pointed out that, for several reasons, 
eq (2) is only approximate. First of all, the integral 
should contain a coefficient 1-^Jb, which is dependent 
on s. Secondly, there should be an obliquity factor 
of the form K(l + cos TPR) [Longhurst, 1957] to 
take account of the fact that the elements of the 
wavefront do not radiate equally in all directions. 
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For the type of applications for which this work is 
intended, however, these details are unimportant 
since the angle a in figure 1 will be relatively small. 
One further objection to the present treatment might 
be that the phase of the received signal is not given 
correctly. This again is not important since it is 
only changes in the phase that are required. How- 
ever, the correct phase could be obtained by per- 
mitting each element of the wavefront to reradiate 

with a phase lead of — relative to the exciting signal. 

These refinements are included in KirchofTs formula- 
tion of Huygen's principle [Longhurst, 1957]. 

Let us now consider an example in which the 
wavefront extends from s to s 3 . Between s 2 
and s 2 the height of the ionosphere is lowered so 
that in this region the phase of the signal is advanced 
by an amount, p [Wait, 1961], but between s Q and 
Si, and s 2 and s 3 , it is unchanged. For brevity let 

E sin ut—Y (0+6) =A 
i?ocosf"cof~ (a+b)~\=B 



and let 



/ abX 

\2(a+T) 



v. 



(6) 



Then using eq (1) and the definitions 

f\os(^)dv=C M> 



(7) 



for the two Fresnel integrals, eq (2) can be written 
in the form, 

E t =A { C{v,)-C{v ) +cos <p[C(v 2 ) - Civ,) ] 

+sin <p[S(v 2 )-S(vO]+C(v 3 ) -C(v 2 ) } 

-BiSW-SivoJ+cosdSivJ-Sfa)] 

+sin pICW-^OI+SW-SW} (8) 



or 
Thus, 

where 



E,=AD-BF. 



E t =^/D 2 +F 2 ■ E sin |~^-y (a+b)-B I (9) 



e=t&w- l FjD 



The change of phase Ad in the received signal is thus 



7T . 7T 

A6=6—- since 6= - when cp=0. 
4 4 



We can now consider some special cases. First of 



all, let <p = Q, which represents the undisturbed case. 
Then, 

D=[-C(v ) + C(v 3 )], 

F=[-S(vo)+S(v z )]. 

To include the whole wavefront, let s and s 3 become 
very large. Then, 



%- 



-oo and # 3 — >-foo . 



Thus, D=l, F=l, and 0=|. Since 

This shows that the total field produced by the 
unperturbed wavefront lags the field due to the 
portion of the wavefront at Q by j and has a value 

V2^ . 

Secondly, consider the case in which v = Vi = — oo 
and v 2 = v 3 = + oo . This represents the case in which 
the whole wavefront is advanced in phase by an 
amount <p. 

Then, 

Z?=cos <£>+sin <p, 

F=cos ^? — sin if. 
Thus, 



?=-—<£>, and A0= — cp, 



wliich shows that the phase of the received signal 

was advanced by an amount tp } as might be expected. 

Thirdly, we consider the more general case in which 



Vq=— oo, V 3 =+co , 



Then, 

D={l-[C(v 2 )-Ctv 1 )](l-cos <p) 

+ [S(v 2 )-S(v l )]sm<p} (10) 
and 

F={l-[Stv 2 )-S(v 1 )](l-cos <p) 

-[CM-CiviXsiiKp]. 

These expressions for D and F, and consequently 
for R and 6, can be evaluated numerically with the 
aid of a table of Fresnel integrals, or graphically by 
means of Cornu's spiral, as shown in figure 2. Con- 
sider a disturbance which extends from .v, to s 2 , where, 
as in figure 1, Si and s 2 are the distances from the 
edges of the disturbance to the line TR. In this case, 
however, it is assumed that the entire disturbance is 
on one side, only, of the line TR. By means of (6), 
the equivalent values v x and v 2) of Si and s 2 are ob- 
tained and marked on the Cornu spiral in figure 2. 
The contribution to the received signal from v= — oo 



711 




Figueb 2. Cornu spiral. 
The received signal is represented in phase and amplitude 1 by JoJ when the 
wavefront is undisturbed. Theeflecl of a phase advance <p on the wavefront 
between^] and V2 is to rotate v% to /':>'. The received signal is now given by JoJ' 
which is the vector sum of the unchanged components i/o»i-H>2 J=Jok, and the 
altered component r^V. 




30°-y- A = 20km 



0.8 1.2 1.6 2.0 2.4 2.8 3.2 3.6 
V = HALF WIDTH OF DEPRESSION 



4.0 4.2 



Figure 3. The effects of a uniform disturbance (for which 
<p='3()°) centered on the path. 
The actual width, 2s, is obtained from eq (7); the numerical values shown refer 
only to the specific case illustrated on the right. 



> Z> 




Ld >- 
CO < 

< -I 



2.0 3 4.0 

V = DISTANCE OFF PATH 

Figure 4. The effects of a uniform disturbance (for which 

^ = 30°) on the received signal as it moves across the path. 

The disturbance is of fixed width. v 2 -v\=2 which in the case illustrated is 316 

km. The actual distances given refer only to the specific case illustrated on the 

right. 



to v=Vi is given by the chord J Q Vi. The contribution 
from the wavefront between Vi and v 2 is given by ViV 2 . 
The effect of the phase advance, <p, along the region 
of the wavefront between v } and v 2 is obtained by 
rotating v x v 2 by the amount, <£>, to the position, v x v' 2 . 
(This follows from (8).) Since the remainder of the 
wavefront between v=v 2 and +°° is given by v 2 J, 
this is added vectorially to Jqv' 2j and J is moved 
to J' . Then the final resultant is given by J J' , 
both in amplitude and phase. In the absence of the 
phase change between V\ and v 2 , the resultant signal 
is given by J J. Thus, the effect of the phase per- 
turbation on amplitude is \J Q J'\/\J J\, and the re- 
sultant phase change is the angle Ad between J J 
and J J'. 

The results of some calculations made in this wa\ r 
are shown in figures 3 and 4. The first of these 
refers to a phase advance which is symmetrically 
situated about the center of the path. The width of 
the wavefront (v = v 1 = —v 2 ) suffering the phase ad- 
vance is the variable. The calculations refer to the 
case where a = 6 = 5000 km, X=20 km and <p = 30°. 
It will be seen that it is not until the width of the 
disturbance (v 2 — v{) approaches about 1.5 that the 
full effect is observed at the receiver. As the width 
increases further, the phase and amplitude of the 
signal at the receiver oscillate about the value which 
it would have if the phase advance extended from 

V= —o° tO +oo . 

Figure 4 shows another example. In this case, a 
phase advance (<p=30°) of fixed width, v 2 — 1\ = 2, is 
considered, but its center is moved across the path. 
The figure shows that as the center of the disturbance 
moves away from the path, the initial perturbations 
decrease from the value given in figure 3. At suffi- 
cient distances from the path, the phase and ampli- 
tude merely oscillate about the values which would 
be observed in the absence of the ionospheric per- 
turbation. The somewhat peculiar shape of the 
oscillation is due to the fact that the Fresnel integrals 
are not periodic for small values of v. 

Figures 3 and 4 have been calculated for the case 
in which the disturbance is equidistant from the ends 
of the path so that a = b. Other cases in which 
ay^b are easily considered by merely decreasing the 
value of s equivalent to a given value of v, in accord- 
ance with (6) . Thus the width of a disturbed region 
of the wavefront necessary to produce a given effect 
is reduced as the disturbance approaches either the 
transmitter or receiver. Furthermore, for off-path 
disturbances, the distance from the path at which a 
particular effect is produced is reduced as the dis- 
turbance moves towards either end of the path. 
Equation (6) also shows that a decrease in wave- 
length has the same effect as moving the disturbance 
towards one end of the path. 

2.1. Effects of Circular Depressions 

In the preceding sections, the profile of the dis- 
turbance in height of the ionosphere has been speci- 
fied by the phase change, <p that would be suffered 
by a ray of infinitesimal width passing through it 
[Wait, 1961]. The effect of varving the width of 
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the disturbed region has then been discussed as 
though the width (s 2 —Si) or (v 2 —Vi) } and <p were 
independent. 

As a first attempt to investigate a more realistic 
case in which (s 2 — $i) and <p are not independent, 
consider a depression which is a square with a side 
length of d when viewed in plan. Then, since p 
increases as the length of path through the dis- 
turbance increases, <p tends to be proportional to d. 
However, as the observed phase change, 6, also tends 
to increase with d, the overall effect, when d is small, 
is for 6 to increase approximately as d 2 . If the depth 
of the depression also increases with d, the depend- 
ence is even stronger. 

In the case of a circular disturbance, there is a 
further complication since the phase change suffered 
by a ray now depends on the distance from the center 
at which it passes through the disturbance. It is 
a maximum for the ray passing through the center, 
but zero for the tangent ray. Furthermore, in most 
cases, the depth of the depression will be a maximum 
at the center and will smoothly decrease until it is 
zero at the edges. Thus the dependence of phase 
shift on distance from the center may become 
marked, and under these circumstances, it might 
be regarded that the "effective" diameter is much 
less than the actual diameter. This effect will be 
illustrated in the next paragraph. 

Wait [1961] has given expressions which enable 
the phase change produced by an infinitely wide 
depression having a vertical profile of the form 

h=h Q for x^>x 

h=h —Ah cos ( - — \ for x<Cx 

to be obtained. In this, h is the height of the un- 
disturbed ionosphere, Ah is the maximum depression, 
x is the distance from the center of the depression, 
and 2x is the width of the depression, both measured 
in the direction of propagation. Wait has calculated 
that such a depression, for which #0=400 km, 
Ah =20 km, h =H0 km, and \=20 km, would ad- 
vance the phase of an infinitesimally narrow ray by 
+26°. Consider now a circular depression whose 
central cross section is of the same cosinusoidal form 
as that just discussed, i.e., of diameter 2x and 
depth Ah. A chord drawn horizontally through this 
depression at a distance x from the center will now 
have a height profile given by 



h'=h —Ah' cos 



( **' \ 

V4V4-Z 2 / 



where x' is the distance from the center of the chord, 
2^x? — x l is the length of the chord and Ah f = Ah cos 

T X 

2'x 

To a first approximation IV is a cosinusoidal 
function of x'. Thus the formulas given by Wait 
[1961] can be used for obtaining the phase delay 



through the disturbance, after calculating the ap- 
propriate values of Ah' and -yjx%— x 2 , for each value 
of x. The results of calculating in this way the 
phase advance for a narrow ray passing through the 
circular disturbance discussed above are given in 
table 1. 



Table 1. Phase change through an 800 km diameter, 20 km 
deep cosinusoidal depression at varying distances from its 
center 



Distance from center. 
Phase advance 



400 
C 



320 

4.4 C 



ii. r 



100 
17.' 



km 
26.2° 



By averaging these phase advances over the width 
of this disturbance, it is found that the same average 
phase advance would be produced by a depression 
only 320 km wide, if its cross section along the 
direction of propagation is everywhere the same as 
the cross section through the center of the circular 
depression. 

The actual effects on the received signal produced 
by a nonuniform disturbance, such as the one dis- 
cussed above, can be determined by using, instead of 
(10), the following expressions. 



-{ 



1 + 0(1^-0^)+^ co S ^[C(v n+1 )-C(v n )} 



f^ sin^' S(v n+1 )-S(:i\)}\ 

71=2 J 



+S 



H 



F=l l+S(v l )-S{v m )+'£; cos^SK+iJ-Sfa)] 



-Ssin^„[C(» B+1 )-C(« 



■„)\y- 



(id 



In these expressions, the disturbed portion of the 
wavefront is divided into m small portions, and </>„, 
is the phase delay produced by the section between 

v n and v n +i. 



2.2. Effects Due to the Curvature of the Earth 

In calculating the phase anomaly caused by a 
cosinusoidal depression in the ionosphere, Wait 
[1961] included the effects of the earth's curvature 
on phase velocity. In the present case, however, 
there is a further geometrical effect which follows 
from the fact that in figure 1 the triangle TPR is 
spherical. If the side lengths are expressed in 

angular measure and denoted by a, 6, etc., then it 
is easily shown, using spherical trigonometry, that 
(1) now becomes 



* A 9 (tan a -ft an b) 

= S A 

2 tan a -tan b 



(12) 
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On substituting 



a / tana- tan b X 
\ 2(tana+tan6) 



(13) 



we again find that ~ir = ~j~ := ~w~ an °l t ne remainder 

of the formulas are unchanged. The consequence of 
this is that s=s • r (where r is the radius of the earth) 
must be increased to obtain a particular value of v. 

Indeed, as a and b both approach 10,000 km( -=- y 

the value of s required for a finite v approaches 
infinity. Under these conditions, the size of the 
disturbance necessary to produce a given effect on 
the received signal becomes very large. It is worth 
noting, however, that a or b separately may approach 
10,000 km since (13) reduces to 



s=- 



X tan b 



V2(l 



+tan a /tan b) 



or 



4 



X tan b 

s=r\ — - v 



when a->-, and similarly for 6->~. Thus when 

either the receiver or transmitter is about 10,000 km 
from the ionospheric depressions, finite results will 
be observed. 

For paths of moderate lengths, the effect of the 
curvature of the earth is quite small. For example, 
consider the case in which a = b = 5000 km, A = 20 
km, and r=6400 km. Then for a flat earth (6) 
shows that s = 158v km, while (13), for the curved 
earth, gives s=180fl km. The differences become 
greater, of course, as a and b increase. 

3. Discussion 

The main purpose of this paper is to enable 
estimates to be made of the effects of a particular 
disturbance as its position relative to the radio 
path is changed. Such estimates are best given 
by comparison with the effects which would be 
observed in the hypothetical case of the received 
signal consisting of an infinitely narrow ray passing 
through the center of the disturbance. The total 
signal at the receiver is the vector sum of two 
portions: that portion which is not altered by the 
disturbance, and that portion which is altered. 
In figure 2 , the vector sum J K of J v t and v 2 J 
comprises the undisturbed portion of the field, 
while v x V2 is the portion which is modified by rotation 
of v 2 to v' 2 through the angle <p. As a general prin- 
ciple, then, from an examination of the Cornu 
spiral in figure 2, it would seem that for the phase 
of the vector sum of the whole wave to differ from 



that due to the wholly undisturbed wave by an 
amount approach ing <p, the len gth ViV 2 =KJ' must 
be greater than J Vx J rv 2 J=J K. This means, in 
essence, that a disturbance centered on the path 
must extend from ^ = 0.7 to #2>+0.7, approxi- 
mately, while the nearest edge of an off-path dis- 
turbance must be close enough to the path that V\ 
be less than unity. These figures are, of course, 
rough approximations and are given only as a 
guide; in the case of off-path disturbances, in partic- 
ular, estimates should be made from the Cornu 
spiral for each special case. 

It should be noted also that there is, in general, 
no simple proportionality between the phase change 
<p of the disturbance, and the phase change Ad 
of the total signal. Examination of figure 2 shows, 
in fact, that the maximum possible phase change, 
Ad', in the received signal is approximately Ad' = 

1 =UL and occurs when ^~J 2 is at right angles to 

J qK. 



tan 



JoK. 



It is perhaps also worthy of mention that if the 
ionospheric disturbance introduces attenuation as 
well as a change of phase, this can be readily taken 
into account. Graphically, as in figure 2, for 
example, the attenuation can be included by simply 
reducing the length of the disturbed vector ViV 2 at 
the same time as it is rotated. 

4. Conclusions 

The effects of localized ionospheric depressions 
on the phase and amplitude of a VLF signal have 
been discussed in terms of elementary Fresnel 
diffraction theory. This has been made possible 
by two main simplifications: (1) by assuming that 
only one waveguide mode is present and that the 
only effect of the depression in the ionosphere is to 
locally alter the phase velocity of the waveguide 
mode, and (2) that the ionospheric depression is 
sufficiently small that the change in phase velocity 
caused by it can be represented as a "lumped" 
phase change in the wavefront where it passes 
through the disturbance. 

It has been shown that: 

(1) Even when the center of the disturbance 
lies along the path, it must be of finite dimensions 
before the phase and amplitude of the received 
signal are affected by it. A small circular dis- 
turbance produces an effect which initially increases 
about as the square of its diameter. 

(2) When the center of the disturbance is not on 
the path, its near edge must pass within a certain 
distance of the path to produce significant effects. 
Otherwise, as the distance of the disturbance from 
the path increases, both the amplitude and phase 
oscillate about their undisturbed values, and phase 
retardations are as likely to be observed as phase 
advances. Indeed, when the center of the dis- 
turbance is at certain distances away from the 
path, the phase of the received signal may show no 
change at all. Under these circumstances, however, 
the amplitude of the received signal will be changed. 
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(3) A disturbance produces the smallest effects 
when it is equidistant from the ends of the path. 
As it is moved nearer one end or the other, its 
effects on the received signal increase. An off-path 
disturbance must, however, move nearer the path 
as it moves towards either end of the path if it is 
to maintain the same effect. It is also easily 
shown that the effects of a given on-path disturb- 
ance are greater at higher frequencies, but again, 
off-path disturbances must be closer to the path 
than at lower frequencies to produce the same effect. 

(4) A geometrical effect due to the curvature of 
of the earth causes the critical dimensions (diameter 
and distance off-path) to be increased above those 
which would be calculated on a flat earth basis to 
produce a change of given magnitude. In particular, 
a disturbance at the center of the path, which is 
about 20,000 km long, will only produce very small 
effects. 

Most of these results seem to be in agreement 
with those obtained independently by Wait [1964 
a, b]. 
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